Two-dimensional double-quantum spectra reveal collective resonances in an atomic 

vapor 
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We report the observation of double- quantum coherence signals in a gas of potassium atoms at 
twice the frequency of the one-quantum coherences. Since a single atom does not have a state at 
the corresponding energy, this observation must be attributed to a collective resonance involving 
multiple atoms. These resonances are induced by weak inter-atomic dipole-dipole interactions, which 
means that the atoms cannot be treated in isolation, even at a low density of 10 12 cm -3 . 



Resonant interactions occur in a wide range of physical 
and chemical systems, from nuclear spin-spin coupling [l[ 
to excitonic effects in light harvesting antennae and re- 
action centers 01. Two-dimensional Fourier transform 
(2DFT) methods excel in isolating and characterizing 
these interactions. Multidimensional Fourier transform 
techniques were originally developed in nuclear magnetic 
resonance, where they proved very powerful in disentan- 
gling congested spectra and providing structural infor- 
mation of complex molecules by measuring the strength 
of inter- nuclear couplings [1]. Thanks to recent devel- 
opments in ultrashort pulse laser technology, these tech- 
niques have been extended into the infrared and visible 
regions of the electromagnetic spectrum during the last 
two decades [4j,|5[. When visible light is used, the method 
is often called 2D electronic spectroscopy because elec- 
tronic transitions are probed. Measurements have pro- 
vided evidence that quantum coherence plays a role in 
photosynthesis [2|, l3[ and given insight into many-body 
interactions in semiconductor nanostructures [6]. The 
ability of 2DFT signals to dissect the quantum pathways 
of matter contributing to the nonlinear response makes 
it a powerful spectroscopic tool. 

Quantum pathways that include a double- quantum co- 
herence can be observed in 2DFT spectroscopy for the ap- 
propriate time ordering of the excitation pulses. They are 
particularly interesting because theoretically they have 
been shown to be highly sensitive to the presence of 
many-body interactions [7[ and can reveal fine details of 
many-body wavefunctions. Double-quantum coherences 
have been observed in 2DFT spectra using infrared ex- 
citation of molecular vibrations [1, @] , in semiconductor 
quantum wells [Io|, EH and in molecular electronic tran- 
sitions [T2I, [Hj . These signals characterize coherences be- 
tween the ground state and a doubly excited state, which 
is often taken as proof that a doubly excited state exists 
in the level structure of an individual chromophore. How- 
ever, we show that collective resonances resulting from 
inter-atomic interactions can induce double-quantum co- 
herences. It has been suggested that such resonances 
can be induced in non-interacting atoms by entangled 



photons |14|, ll5|- Interaction-induced double-quantum 
coherences may play a role in photo-processes such as 
light harvesting [2|, [3|. Furthermore, quantum informa- 
tion processing schemes have been proposed that exploit 
collective resonances due to long-range dipole-dipole cou- 
pling [l6j. 

Here, we report the observation of double-quantum 
coherences in a gas of potassium atoms. Based on the 
level scheme of a single potassium atom, shown in Fig. 
1(a), and the spectrum of the excitation laser, the naive 
expectation is that double-quantum coherences should 
not exist because there are no atomic levels at twice the 
frequency of the laser, which is tuned to be resonant 
with the one-quantum transitions corresponding to the 
Di and D2 lines. We interpret the observed two-quantum 
resonances as being due to collective excitations of two 
atoms, and use a Hilbert space transformation to provide 
a heuristic explanation of how interactions result in an 
effective level at twice the one-quantum energy. For sim- 
plicity, we first discuss two- level atoms, as the essential 
physics is not changed because we excite both Di and D2 
transitions. A more rigorous theoretical treatment using 
an excitonic picture and incorporating dipole-dipole in- 
teractions in a gas of atoms combined with molecular dy- 
namics simulations reproduces the experimental results. 

The quantum pathways that include a double- 
quantum coherence are distinguished by the combina- 
tion of the wave-vector and time-ordering selection of 
the pulses. As shown in Figure 1(b), the three incident 
pulses, with wavevectors k^, k# and kc, focus and inter- 
sect at the sample such that they are on three corners of 
a square that is perpendicular to the propagation direc- 
tion. Their interaction produces a signal in the direction 
ks = — k^ + k# + kc, which corresponds to the 4th 
corner of the square after the sample. Double-quantum 
coherences can contribute to the signal only if the pulse 
with the conjugated wavevector, k^, arrives last. This 
excitation sequence is shown schematically in Fig. 1(c) 
and (d) for a simplified 3-level system consisting of a 
ground state, |0), singly excited state, |1), and double 
excited state, |2). The first pulse, kg, puts the system in 
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FIG. 1. (color online) (a) Levels of an isolated potassium 
atom; (b) Geometry of incident beams and signal beam; (c) 
Pulses and oscillating coherences during different time peri- 
ods; (d) Transitions of a ladder-level system as a function of 
time showing transitions driven by laser pulses and result- 
ing coherences; (e) Double sided Feynman diagrams for the 
two quantum pathways that contribute to Sm for a ladder 
level scheme; (f) Hilbert space transformation between two 
independent two-level atoms and a 4-level system including a 
doubly excited state. 



a quantum mechanical superposition between states |0) 
and |1), i.e., in a single quantum coherence. The second 
pulse converts the single quantum coherence to a dou- 
ble quantum coherence between |0) and |2). The third 
pulse converts the double quantum coherence back to a 
single quantum coherence that radiates. The final single- 
quantum coherence can be between states |0) and |1) or 
between states |1) and |2), both of these possibilities are 
shown in Fig 1(d). These pathways can be described us- 
ing the two double-side Feynman diagrams for the atomic 
density matrix shown in Fig. 1(e). 

Experimental two-dimensional single- and double- 
quantum spectra are acquired using the apparatus de- 
scribed in Ref. [13. A mode-locked Tiisapphire laser gen- 
erates ~ 200 fs pulses that are input to an ultrastable 
platform of nested and phase stabilized interferometers 



to generate 4 identical pulses arranged in a box geometry. 
Three of the beams are the excitation beams, k^, kg and 
kc, while the fourth is designated the tracer and propa- 
gates in the same direction as the signal beam, kc- The 
tracer is used to generate a reference pulse that is routed 
around the sample and interfered with the emitted sig- 
nal beam to produce inter ferograms. The full phase and 
amplitude information about the signal can be extracted 
from the inter ferograms. The tracer itself is blocked dur- 
ing data acquisition. It is used to determine the overall 
phase of the signal for decomposition into real and imag- 
inary parts [18] . 

The potassium vapor is held in a 350 micron thick 
transmission cell. The cell body is made of titanium 
with two sapphire windows that are diffusion bonded to 
the titanium. The vapor temperature is controlled by 
a heater attached to the cell. For the measurement re- 
ported here, the cell temperature was 130 C. The trans- 
mitted spectrum of the strongly attenuated laser was 
used to estimate the absorbance. The transmitted inten- 
sity is I t = Iie~ a \ where 1{ is the incident intensity, a is 
the absorption coefficient and I is the cell thickness. At 
this temperature and density, the resonance broadened 
linewidths are smaller than the Doppler width and spec- 
trometer resolution, thus a 1550 torr Argon buffer gas 
was introduced into the cell to induce collision broaden- 
ing and reduce the peak absorbance al to 0.053 for the 
D 2 line, which is stronger than the D\ line. This low 
absorbance rules out optical density effects as explaining 
the observations. 

Figure 2 shows the real part of both single- and double- 
quantum spectra for a potassium number density of 
3.5 x 10 12 cm -3 . As observed in previous measurements 
[HI, the single-quantum spectrum shows peaks corre- 
sponding to the Di and D 2 lines, as well as off-diagonal 
peaks due to coupling between them via the ground- 
state bleach and Raman-like coherences. Surprisingly, 
the double-quantum spectrum also shows clear peaks 
even though there are no atomic states at these ener- 
gies. The observed linewidths correspond to a dephasing 
time slightly longer 5 ps for the double-quantum coher- 
ences. The Di and D 2 lines correspond to transitions 
from the 4 2 S 1 / 2 ground state to the 4 2 P 1 / 2 and 4 2 P 3 / 2 
states, with transition frequencies of 389.29 and 391.02 
THz, respectively. The higher lying states that are clos- 
est to twice the Di j2 energies are the 5P and AD states, 
which are at frequencies 740.81 and 821.36 THz, both of 
which are well outside the spectral range of shown in Fig. 
2(a) and well outside the laser bandwidth of 3 THz. The 
fact that the observed resonances are at exactly twice 
the frequencies of the D\ and D 2 lines (or their sum fre- 
quency) indicate that the double-quantum resonances are 
due to the combined response of two atoms, rather than 
the level structure of a single atom and that the coupling 
is weak (otherwise the resonance would be shifted). 

The observations can be explained using a simple pic- 
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FIG. 2. (color) Experimental (a) double-quantum and (b) 
single-quantum real spectra. 



ture obtained by combining the Hamiltonians of two in- 
dividual atoms into a Hamiltonian describing both si- 
multaneously, as shown in Fig. 1(f), and working in 
the joint Hilbert space of both atoms. For two isolated 
atoms, each described as a two level system, there are 
no doubly-excited levels, so a double-quantum coherence 
is impossible. However, the joint space of the two two- 
level atoms has 4-levels. These consist of a ground state, 
where both atoms are in the individual ground states, 
two singly excited states where one of the two atoms is 
excited, and a doubly excited state where both of them 
are excited. In this four-level system, double-quantum 
coherences can occur between the ground state and the 
doubly excited state. However, this change of description 



should not affect the underlying physics, which seems 
contrary to the fact that the nonlinear susceptibility cal- 
culated for the four-level system has terms that are not 
present for two isolated two- level systems. Careful analy- 
sis shows that all of these "new" terms for the four-level 
system interfere destructively and exactly cancel if the 
ground state to singly-excited state transitions are iden- 
tical to the singly-excited to doubly-excited state tran- 
sitions. For the double- quantum coherences, this means 
that the two pathways shown in Fig. 1(e) cancel one an- 
other. Thus describing the system in the joint Hilbert 
space alone does not change the nonlinear susceptibility. 
The picture changes if there are interactions between the 
atoms. In that case, the transitions will not be identi- 
cal and the cancelation will be incomplete, resulting in 
double-quantum resonances. The joint Hilbert space de- 
scription is useful in explaining why interactions result 
in double- quantum resonances. This picture only consid- 
ers the interaction of a pair of atoms; a more complete 
theoretical treatment needs to consider an interacting en- 
semble. 

To go beyond this simple explanation, we performed 
simulations using the exciton formalism. Excitons are 
delocalized collective electronic excitations of assemblies 
of coupled chromophores including coherences that re- 
flect entanglement. Excitons with extended coherence 
size were first discussed in crystals 20|, El|. In real crys- 
tals, disorder and dynamic fluctuations limit the coher- 
ence size. An atomic vapor is an ideal model system for 
studying excitonic effects in the presence of disorder [22 1. 
In atomic vapors, n-exciton states correspond to states 
where n atoms are excited and all others are in the ground 
state. The atoms interact due to dipole-dipole coupling. 

The single- and double-quantum spectra are simulated 
based on the Nonlinear Exciton Equations (NEE) 
The effective Frenkel exciton Hamiltonian, which de- 
scribes the electronic excitation of interacting atoms, is 

H(t) = Y J ^m{t)Bl n B m + 

m 

J2 (Jmn(t)BlB n + K mn (t)BlBlB m B n ) 

(i) 

where B^ and B m are the exciton creation and annihila- 
tion operators on state m, m is a multi-index consisting 
of the atom number and the intra-atom exciton state(e.g. 
states connected to D\ or D2 line). Q m is the excitation 
energy of exciton m, J mn indicates the one-exciton cou- 
pling strength between excitons m and n and K mn is 
the two-exciton coupling parameter. The dipole-dipole 
coupling model for J mn can be calculated from 



Jmn (t) 



1 



3(/i m • Rmn(t))(P>n ' Rmn(t)) 



1 L mn\ b ) 



(2) 



4 



where j2 m ,n are the transition dipole moments and R mn 
is the distance between atoms m and n, e is the dielectric 
constant. The Hamiltonian for the coupling of the atoms 
to the optical field is 

H int (t) = -E(r,t)-V(t) 

= ~E{v,t) ^itrn(t)(Bl(t) + B m (t)). ^ 

m 

All elements of the Hamiltonian are time dependent 
through the changing configuration of the potassium 
atoms. 

The simulation of the double-quantum coherence spec- 
tra is based on a direct nonlinear exciton propagation 
method in an exciton picture [23j. The simulation of 
the 2DFT spectra requires three steps: first, the trajec- 
tories of an ensemble of atoms are simulated; second, 
the trajectories of the atoms are used to calculate the 
time-dependent interatomic couplings for atoms; third, 
the 2DFT spectra are generated using the calculated in- 
teratomic couplings. Dynamic trajectories of 20 potas- 
sium atoms in a simulation unit cell were generated us- 
ing the molecular dynamics package GROMACS. The 
ground state potential energy surface of K 2 is used as 
the coupling function between atoms. Runs up to 2 ns 
were performed with a resolution of 0.01 ps. Then the 
trajectories of atoms are used to calculate the interatomic 
distances between atom pairs and the coupling parameter 
J mn for each snapshot of the trajectory using Eq. 2. The 
intra-atomic elements of K mn are set to large numbers to 
ensure that each atom cannot be excited twice, while the 
interatomic elements of K mn are set to zero since the in- 
teratomic two-exciton couplings are ignored. As the last 
step, 2D spectra are generated using the SPECTRON 
package [23[. The double-quantum coherence spectrum, 
S///(r,T,t),is 

S in (r, T, t) = 2 (^) 4 &n ( r + T + 1 + **) 

t s ri4 

xRnl 11 (r + T + t + t s , T + t + t a , t + ts,t a ). (4) 

The signal is averaged over several start times, t 3 , along 
the trajectory with 

dsG n4 , m4 (T 4 , s)K m47ni (s) 

7/7,47/7,1 3 

^ImM^n^is;^) (5) 

for the one-exciton Green function G n2jni (r2,ri) = 
( ff |B n2 t/(T 2 ,Ti)Bt 2 | ff ) where U(t 2 , ri) is the time- 
evolution operator connected to H(t) . The excitonic 
functions are computed by direct integration of the 
Schrodinger equation ihd t \^^ (t; T3)) = H (t)^ 1 )^; T3)) 
and ihd t \ilJ i2) (t;T3)) = H (£)|^ (2) (t; r 3 )> with the initial 



conditions ipmiit = t\\t\) = /i mi (ri) • E(t\) and 

^m!,m 2 (* = ^2;r 2 ;ri) = 

-^{Mm 2 (r 2 ) 'E(r 2 )^l(t = r 2 ;ri) 
H-Mmx^) • ^(r 2 )^(t = r 2 ;n)}. 

The double-quantum spectrum shown in Fig. 3(a) is 
obtained by Fourier transforming Siu(r,T,t) with re- 
spect to r and t. The single-quantum spectrum (Fig. 
3(b)) was calculated similarly as for the double-quantum 
spectrum. The model parameters have been adjusted to 
match the experiment; however our goal is not a detailed 
fit to the experimental results, but rather to simply show 
that double-quantum resonances do appear. The simula- 
tions do show double-quantum coherences, which do not 
appear if the dipole-dipole coupling is turned off. Fur- 
thermore, the simulated lineshapes are in good agreement 
with the experiment, for example both display a disper- 
sive profile, which has been shown to be characteristic of 
many-body interactions in semiconductors [24[. 

Our results demonstrate that, thanks to its high sensi- 
tivity, electronic double-quantum coherence spectroscopy 
can detect even very weak interactions between chro- 
mophores, which makes it a powerful technique, but at 
the same time means that the interpretation of double- 
quantum spectra may not be possible based on the level 
structure of a single chromophore and thus collective ef- 
fects should be considered. Our results are for the regime 
where any spectral shifts are small compared to linewidth 
and laser bandwidth, as compared to observations fo a 
Rydberg blockade that were in the opposite regime [25|. 
Moreover, the long-range nature of the dipole-dipole in- 
teraction means that this effect is present even at low 
densities, where it is often thought that interactions can 
be neglected. Collective coupling to the radiation field 
is also responsible for enhanced spontaneous emission or 
superradiance [26] . We find that the ratio of the strength 
of the double-quantum coherences to the single-quantum 
coherences is constant as the atomic density is varied over 
an order of magnitude. This result is consistent with the 
fact that self-broadening is theoretically found to be in- 
dependent of density [27] due to the r~ 3 scaling of the 
dipole-dipole coupling. 
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FIG. 3. (color) Simulated (a) double- quantum and (b) single- 
quantum real spectra. 
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